Evidence for multiferroicity in TTF-CA organic molecular crystals 
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We show by means of ab-initio calculations that the organic molecular crystal TTF-CA is mul- 
tiferroic: it has an instability to develop spontaneously both ferroelectric and magnetic ordering. 
Ferroelectricity is driven by a Peierls transition of the TTF-CA in its ionic state. Subsequent anti- 
ferromagnetic ordering strongly enhances the opposing electronic contribution to the polarization: it 
is so large that it switches the direction of the total ferroelectric moment. Within an extended Hub- 
bard model we capture the essence of the electronic interactions in TTF-CA, confirm the presence 
of a multiferroic groundstate and clarify how this state develops microscopically. 

PACS numbers; 71.45.Gm, 71.10.Ca, 71.10.-w, 73.21.-b 



Introduction The search for multiferroics - single phase 
materials in which magnetism and ferroelectricity coex- 
ist - has become an important research topic in the last 
few years A large number of new multiferroics have 
been discovered, often guided by the theoretical predic- 
tion on the presence of multiferroicity Almost all 
newly discovered multiferroics are transition metal com- 
pounds where spin, lattice and charge degrees of freedom 
are strongly entangled. Here, we explore a new direction 
and extend this search towards organic systems. Synthe- 
sizing organic ferroelectrics is a rapidly evolving endeavor 
within the field of organics, materials that are known for 
their lightness, flexibility and non-toxicity Finding 
an organic multiferroic can open up a new area of ma- 
terials research with potentially mechanisms for the si- 
multaneous occurence of magnetic and ferroelectric order 
different from those active in standard transition metal 
oxide multiferroics. Here we show with a combination 
of density functional theory (DFT) and model Hamil- 
tonian calculations that the organic molecular crystal 
TTF-CA is ferroelectric. The magnitude of the polar- 
ization strongly depends on the presence of magnetic or- 
der, thereby providing an interdependence of magnetism 
and ferroelectricity. To the best of our knowledge this 
is the first theoretical prediction for multiferroicity in an 
organic material. 

Neutral-Ionic Transition TTF-CA (tetrathialfulvalene- 
p-chloranyl) is a charge transfer (CT) salt with stacks of 
alternating donor (TTF) and acceptor (CA) molecules. 
This material has been studied over the last few decades 
because of its prototypical neutral-ionic transition (NIT) 
At ambient pressure, TTF-CA undergoes a first or- 
der NIT at a critical temperature of ~ 81 K[7|. The crys- 
tal symmetry is lowered from P12i/7il (centrosymmet- 
ric) for the neutral state to Plnl (non-centrosymmetric) 
for the ionic phase (il . This symmetry breaking changes 
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FIG. 1: Schematic arrangment of TTF and CA molecules 
along the a axis in a) the neutral (N) b) ionic (I) structural 
phases with finite spontaneous polarization, where arrows in- 
dicate the displacement of CA towards TTF due to the dimer 
formation, and c) AF magnetic structure in the ionic phase. 



the position of TTF and CA molecules within the plane, 
causing the formation of pairs of long and short bonds 
along the stacking axis a (see Fig. [T|). The formation 
of dimers is reflected in the dielectric response, which 
shows that the NIT in TTF-CA is associated with a 
first-order ferroelectric phase transition [i,[l^. Both the 
non-polar and polar states of TTF-CA are electronically 
charged 11. 1^. 

Ab initio electronic structure We perform density func- 
tional theory calculations [13] using the generalized gra- 
dient approximation to exchange-correlation functional 
according to Perdew-Becke-Erzenhof (PBE) [l3|. The 
electronic structure is computed using the projector aug- 
mented wave method (PAW) [li| , as implemented in the 
Vienna ab-initio Simulation Package (VASP)[il, [13] ■ 
The electronic contribution to the electric polarization 
is calculated by the Berry phase (BP) method developed 
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FIG. 2: Non-spin polarized band structure of Plnl experi- 
mental structure of Ref. (0] along high simmetry directions 
(labelled as: X=(i,0,0), S=(i,0,|), Z=(0,0,|), r=(0,0,0), 
Y=(0,i,0), U=(0,i,i), T=(i,i,i), W=(i,i,0)) and density 
of states (DOS) projected onto TTF and CA molecule and 
calculated using either PBE (left) or HSE (right) functional. 



by King-Smith and Vanderbilt [181 . 

Calculations for the high temperature centrosymmetric 
structure (space group P12i/ nl |8|]) provide a simple pic- 
ture for the electronic structure: the TTF-HOMO (High- 
est Occupied Molecular Orbital) and CA-LUMO (Lowest 
Unoccupied Molecular Orbital) of the four molecules hy- 
bridize, giving rise to two bonding and antibonding bands 
around the Fermi level. This situation in essence corre- 
sponds to an half-filled band system. Electronic charge 
is donated from TTF to the CA, through occupation of 
the bonding bands The bandwidth of the occu- 

pied manifold is Vl^=1.0 eV, the system is metallic and it 
shows a quasi one-dimensional character along the stack- 
ing axis a as is highlighted in the band structure results 
by the dispersion along the F-X direction of the Brillouin 
zone [13, [2^. 

Using the low-symmetry structure of TTF-CA (space 
group Plnl [1]) the system opens a very small gap A ~ 
50 meV, with the formation of dimers leading to a neg- 
ligible change in the bandwidth [21*1 . The corresponding 
PBE band structure is shown in Fig. [51 Semi-local func- 
tionals fail to describe correctly some of the most chem- 
ically significant orbitals of organic molecules, in partic- 
ular the HOMO and LUMO molecular orbitals. This 
defect can partly be traced back to a self-interaction er- 
ror, i. e. the spurious Coulomb interaction of an electron 
with itself [22]. Hybrid functionals, which mix a fraction 
of the Fock exchange with the exchange density func- 
tional have shown to largely overcome these shortcom- 
ings for molecules [2^ and molecular crystals We 
thus expect with hybrid functionals a better physically 
description of molecular crystals, also in the light of their 
cumulative "track record" in organic chemistry Q^. We 



note that the DFT-I-U method [27|, often used to deal 
with strong on-site Coulomb correlations in transition 
metal oxides, is not directly applicable in molecular crys- 
tals. The reason is that in molecules the localized orbitals 
are multi-center rather than single-center, since the basis 
set correspond to orthonormal molecular orbitals instead 
of orthonormal atomic orbitals. Here we used the Heyd- 
Scuseria-Ernzerhof hybrid functional (HSE) [28|] which 
does impressively well for extended solid state systems 



As opposed to PBE results, the Non-Magnetic (NM) 
HSE ground state of the ionic state Plnl already has a 
gap: A=0.4 eV[23|. The corresponding band structure 
obtained with HSE is shown in Fig. [2] We also performed 
spin polarized HSE calculations and find that the sublat- 
tice of TTF and CA molecules stacked along the a axis 
become antiferromagnetically (AF) ordered, with a rnag- 
netization of 0.4 hb per molecule (see Fig. [T] c)) [31j. 
The energy gain of the magnetic state with respect to 
the non-spin polarized state is 80 meV per unit cell [s^ . 
The electronic gap in the magnetic groundstate is A=0.5 
eV. 

We evaluate the ferroelectric polarization P in both the 
NM and AF states. We characterize the adiabatic path 
of ionic displacements from the paraelectric (P12i/nl) to 
the ferroelectric (Plnl) structure by A. The parameter A 
is therefore a measure for the amount of dimerization of 
TTF and CA molecules (see Fig. [T]). In the NM state, we 
find a net polarization of P = 8.0/iC/cTO^, which reduces 
to P = S.diiC/cm^ in the AF state (see Fig. The 
polarization is in the a — c plane [s^. Issjl (see Fig. [3^)), 
i.e. P = PaaL + PqC, with an almost negligible Pc compo- 
nent. It is remarkable that Pa depends strongly on both 
the spin state and A: P^^ is much larger than P^^. The 
decomposition of the total polarization into an ionic and 
electronic part, P = P"^'"^ + piomc^ elucidates this obser- 
vation, see Fig. [3)3). These two components both increase 
in magnitude along the adiabatic path (i.e. as a function 
of A) but are always opposite to each other. Such a par- 
tial cancellation of electronic and ionic polarizations is 
also encountered in some inorganic materials, for exam- 
ple in the proper ferroelectric BaTiOs, and in improper 
multiferroic materials such as HoMn9 0»^(36|. where the 
polarization is magnetically driven. 

Compared to the NM situation, in the AF state P^'^ 
becomes much larger due to superexchange driving the 
formation of a larger dipole moment in the electron cloud. 
As P^'-'^ is opposite to the p^°"'= (not aflFected by the 
magnetic ordering), an increase in the AF state of the 
former first causes a reduction in the magnitude of P, but 
when it becomes larger than p^o"*^^ j-j. ultimately causes 
the magnetic ordering to flip the direction of total P (see 
Fig. [li)) - this is what we find to occur in TTF-CA. 

Model Hamiltonian In order to understand the mech- 
anism leading to structural and ferroelectric transitions, 
we model TTF-CA by an extended Hubbard model on a 
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FIG. 3: (a) HSE Ferroelectric polarization along a and c axis 
{Pa and Pc) for AF and NM states along the path connecting 
the P12i/nl (A=0) and Plnl (A=l) crystal structure, (b) 
and (c) Decomposition of P into electronic {Pde) and ionic 
(Pionic) contributions for the NM (panel b) and AF (panel c) 
states. 



one-dimensional lattice with the Hamiltonian 37, 38, 3^, 
+(5^(-irn, + C/^n,tn,X + 1/2^^2. (1) 



Here, Cj^. and c]^ are annihilation and creation operators 
for electrons at site i with spin a, t denotes the bare hop- 
ping strength and Ui are the distortions along the chain 
of the TTF and CA ions from their equilibrium position. 
The electron-phonon coupling parameter is a, which fa- 
vors Peierls dimerization and 5 is the magnitude of the 
ionic potential arising due to inequivalency of the TTF 
and CA and U denotes the Hubbard repulsion strength. 

The hopping parameter is i = 0.2eV as extracted from 
a fit of the PBE band structure to one dimensional tight- 
binding chain and the on-site difference is (5=0.06 eV 
Til 20j . The Hubbard U is the screened energy cost re- 
quired to double occupy a molecular site and can be eval- 
uated from total energy of molecular ions TTF and CA 
and the polarization energy (ill, [i^ lisllii. lislli^ |47|. 
We calculated the variation of the eigenvalues of HOMO 
of TTF and LUMO of CA and their polarization energy 
of the charged molecules placed inside a cavity of an ho- 
mogeneous dielectric medium with dielectric constant of 
3 using the SS(V)PE model We estimate Uttp'=2.3 
eV. which is in good agreement with the value of Ref. 
and UcA=2.7 eV 0, SI, SI. Note that this pro- 
cedure does not account for the role of hybridization of 
HOMO and LUMO and dipole formation of the differ- 
ent molecular TTF and CA sites; it thus gives an upper 
bound on the value for the strength of electron-electron 
interactions. Moreover close to NIT the screening is en- 
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FIG. 4; (Color online) (a) Phase diagram showing undis- 
torted to dimerized, and NM to AF transitions. The shaded 
region on left (right) indicates the ferroelectric (multiferroic) 
state, (b) lonicity parameter p — uttf — ncA as a function 
of U for different values of a. (c) The polarization computed 
within point charge model, and the magnetic moment as a 
function of U for the same values of a as in (b). 



hanced by the diverging dielectric constant d, de- 
creasing the effective Coulomb interactions. In this re- 
spect, we remark that, by performing accurate GqWo 
calculations ^5 1] for the crystal structure of TTF-CA in 
the Plnl symmetry we obtain an upper-bound for the 
quasiparticle gap of 1 eV S^], pointing to a value of U 
considerably smaller than that estimated above. All this 
suggests the importance of hybridization and screening 
in the mechanism driving the NIT. 



Employing a Hartree-Fock decomposition for the Hub- 
bard term, we solve the above model self-consistently 
for different initial conditions and retain the lowest en- 
ergy solutions. Deviation from the undistorted state, i.e. 
dimerization, is realized via a staggered pattern of distor- 
tions Ui. Therefore we set: Ui = mo(— 1)* and treat uq as 
a variational parameter. In Fig. 0^, we show the ground- 
state phase diagram in the parameter space of a/t and 
U /t. On basis of the DFT calculations we fix the value of 
5 to 0.3i. The phase diagram essentially consists of two 
phase boundaries crossing each other: (i) non-magnetic 
ionic state to AFM ionic state with increasing U , and (ii) 
an undistorted state to a Peierls dimerized state with in- 
creasing a. The ionic -|-dimer state is ferroelectric, shown 
as the shaded region in the phase diagram. The ionic- 
ity p, here defined as the difference between the charge 
density at the TTF and the CA sites, arises due to the 
presence of the on-site energy difference 5. It decreases 
monotonically with increasing U (see Fig. |4}d). The ionic 
state with dimerization of TTF-CA lacks inversion sym- 
metry and is therefore ferroelectric. We plot the polar- 
ization calculated within a point charge model and the 
local magnetic moment as a function of U in Fig. 0};. The 
polarization vanishes in the undistorted state due to the 



4 



inversion symmetry. For a = 0.4 and a = 0.6 polariza- 
tion reduces with increasing U in the non-magnetic state 
and shows a sudden drop at the transition from a non- 
magnetic to AFM state due to a reduction in the dimer 
distortion uq. The magnetic moment shows a behavior 
that correlates with the polarization. This presents an 
interesting case where the magnitude of the polarization 
is controlled by the nature of the magnetic state. We find 
a gap A - 2t in the DOS for S/t = 0.3, U/t = 3, and 
a/t = 0.4, which is fully consistent with DPT. 

In conclusion, we performed first-principles calcula- 
tions to understand the electronic and ferroelectric prop- 
erties of the TTF-CA molecular crystal. We find that 
electronic correlations induce magnetism, which in turn 
control the ferroelectric polarization in the material. The 
polarization is present even in the non-magnetic state; 
however, both the magnitude and direction of the polar- 
ization are very sensitive to the magnetic ordering. Using 
an extended ID Hubbard model, we show that indeed, 
both the ferroelectric and the multiferroic states can co- 
exist in the model and the calculation of the polarization 
within a point charge model qualitatively agrees with the 
results of first-principles calculations. Our theoretical re- 
sults and predictions provide guidance for future exper- 
imental investigations on multiferroicity in the organic 
molecular crystal TTF-CA. 
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